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Abstract: The order parameter profile between coexisting confined and 
plasma regions at the QCD phase transition is constructed. The dimen- 
sionless combination of the surface energy (or) and the correlation length 
(£) is estimated to be 

<rf 3 « 0.8. 
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ORIGINAL PAGE IS '' “** 

OF POOR QUALITY 

The impact of the transition from the quark-gluon plasma to the hadronic 
phase on the baryon density fluctuations and primordial nucleosynthesis in the 
Early Universe is being actively investigated [1-4]. Quantitative description of the 
hadronic bubble formation, their percolation and the squeeze of the plasma droplets 
rely on data originating from quantum chromodynamics. The progress in the lattice 
measurements of the transition temperature and the latent heat is continuous and 
at present a rather weakly first order transition is suggested to take place [5,6]. 

The average size, the rate of production etc. of hadronic bubbles are sensitive 
to a third quantity, the surface energy associated with the interface between coex- 
isting confined and deconfined regions at the transition temperature. Very recently 
attempts have been made for numerical determination of this quantity directly in 
QCD [7,8]. The problem consists of calculating the excess free energy of the gluonic 
medium at T c resulting from the imposition of two parallel, infinitely distant walls 
of temperature T = oo and T = 0, respectively. The exploratory Monte Carlo 
study [8] does not correspond yet to this ideal setup. 

In this note we propose as an approximate realisation of the above strategy the 
evaluation of the surface energy in the framework of the mean field approximation 
to the order parameter (Polyakov-loop) theory. The strategy is to construct and 
quantitatively characterize the so called order parameter profile [9]. This approxi- 
mation is widely used in various fields of statistical physics e.g. nucleation theory, 
wetting etc. The explicit form of the effective theory of Polyakov-loops we shall use 
is valid at strong gauge coupling only, still we believe that predictions on combina- 
tions of observables which are independent of the lattice spacing might be of more 
general validity as an order of magnitude estimate. (This expectation is based on 
the experience gained with strong coupling spectrum calculations [10]). 

The derivation of the effective lattice theory of Polyakov-loops should be con- 
sidered as a constrained integration in the partition function of the full, finite tem- 
perature gauge system. Were the integrations performed exactly, the ’’integrated 
out” contributions to the free energy kept, then no thermodynamic information 
should be lost. Furthermore the pieces independent of the Polyakov loop variables 
cannot produce the necessary singularities, therefore the latent heat, the order pa- 
rameter discontinuity etc. should not depend on them. 

To the present date, the integration has been performed only in the strong 
coupling regime [11,12]: 

S e ff = J e ff Ux Tr U+ +e + c.c.), (1) 

*,e 


where 



( 2 ) 


(g 2 is the gauge coupling 1 , T the temperature and a the lattice constant in 
the temperature direction). We consider the pure 517(3) gluon theory. Closer to 
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the continuum limit (g 3 —* 0,a — » 0) more complicated terms do appear in (1), 
with further J-type couplings depending on g 2 and T nontrivially. For individual 
physical quantities there is no hope to see the continuum behavior. The best, the 
predictions for a-independent ratios can be taken seriously. 

Let us denote the direction orthogonal to the boundary walls as the x-axis. The 
physical situation suggests to choose an x-dependent mean field m(x) to minimize 
the approximate free energy per unit area given by the expression: 

Fm.f. = / iV «p{ JM(x)(Tr V, + Tr V+)} + J £ m(x)M(z), (3) 

* " X 


where M(x) = 4 m(x) + m(x + 1) + m(x — 1) and dU is the Haar-measure of the 
SU’(S)- invariant integration. One has to solve the selfconsistent equation: 


6m(x) 


= 0 


(4) 


subject to the boundary conditions 


m(x — + — oo) = 0, m(x — ♦ oo) = mo. (5) 

The parameter mo is selected by the requirement of the existence of a soliton solu- 
tion to eq.(4), smoothly extrapolating between the ordered and disordered phases 
(see below). Eq.(4) is understood as a local relationship between M and m, or 
introducing the notation 


m(x + 1) = m(x) + i(x), 


M{x) = 6m(x) + i(x) — i(x — 1) = 6m(x) + I(x) (6) 


a relation between I and m. One realizes that J(x) approximates cPm/dx 2 , there- 
fore eq.(4) is actually the discrete version of a second order nonlinear differential 
equation. (For instance, in the 3-state Potts model this equation can be written 
in fully explicit form. With this identification the smoothness of the interpolation 
between the prescribed boundary values is naturally expressed by the boundary 
conditions: 


dm 

dx 


-» 0 , 


(Pm 

dx 2 


— ► 0 , 


if |x| — > oo. 


(7a) 


On the lattice these conditions are equivalent to: 


ati»o 

Jo 


I(m)dm = 0, 


J(x) — ► 0, if |x| — ► oo. 


(76) 
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What are the consequences of (7b) on the parameters mo and J't This question 
is answered starting from eq.(4), whose explicit form is: 


1 d 

2 JdM 


In J dUexp{JM(Tr U + Tr U + )} = m. 


( 8 ) 


Integration of eq.(8) with respect to M on the interval (0,M(mo)) results in the 
following equality: 


In 



exp{JM(mo)(Tr U + Tr U + )} — Jm 0 M(m 0 )— 


pm 0 

— mol(mo) — / dm J(m) = 0. (9) 

Jo 

In view of (7b) the last two terms on the j hand side of (9) vanish. In this 
way eqs.(9),(8) are just the equations whi ctermine the transition point J c and 
the order parameter discontinuity (m c ) in the conventional treatment of first order 
transitions. 

The meaning of this result is that eq.(8) has a smoothly interpolating solution 
only at J = J c with mo = m c . Clearly, the physical statement of phase coexistence 
is equivalent to the existence of a solitonic order parameter profile. 

The numerical construction of the kink solution starts by solving eqs.(7b),(8), 
(9) with m(x) = m c , which reproduces the well-known J c = 0.134270 and m c = 
1.45786 data [9]. Using these values in (8), m(x+l) is calculated from the knowledge 
of m(x — 1) and m(x). The ’’initial” values were chosen to be m(0) = 0,m(l) « 
o(lO -10 ). Each step of calculating the soliton in a new point involves two integrals 
to be done numerically. The accumulation of numerical errors might spoil the 
soliton.(The inaccurately calculated m(x) function would reach m 0 = m c with 
nonzero derivative, making the kink unstable). Therefore J c , m c and the integrals 
require careful evaluation. (In case of the 3-state Potts model, the homogeneous 
mean field problem can be solved analytically and the sensitivity to deviations from 
J c and m c be studied in a systematic way). The kink constructed numerically is 
displayed in Fig.l. 

The soliton can be characterized by two physical quantity. Its half-width pro- 
vides an estimate of the order parameter correlation length at the transition. From 
Fig.l one reads off £ « 2 in lattice units. Substituting m(x) into eq.(3) the surface 
energy in lattice units is 

a = 0.106. (10) 

A combination which is independent of the actual lattice spacing can be formed as 


<r£ 3 « 0.848. 


( 11 ) 


This is the main message of our note, strictly speaking valid in the strong coupling 
regime. 
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However we cannot resist to conjecture a value of <r in physical units extending 
the validity of (11) to weaker couplings. For this we use data from the most recent 
Monte Carlo investigations [5,6]. In Ref. [5] a detailed analysis of the order param- 
eter correlation length was given. For ’’temporal size” N t = 4 two different spatial 
sizes led to £ « 9 and £ « 13, respectively at the transition point. Assuming that 
this system (fl c = 5.69) already obeys continuum physics one deduces 

(phyT c » 2.75 ± 0.5. (12) 

Combining this with our result (11) one gets 

<r P h v , * (0.04 ± 0.02)T*. (13) 

Our result for <r p hy/T* is almost two orders of magnitude smaller, than that of 
Ref. [8]. We consider this shift towards smaller coefficients physically reasonable as 
one deals with a weakly first order transition (in continuous transitions <r(T e ) = 0). 
Assuming T e = lOOAfeV" the range of <r p hyt from (13) practically coincides with 
the values considered by Fuller et al. [3] for astrophysical applications. For T c = 
150-MeV the values from (13) are at the higher end of the range considered in [3]. 
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Figure Caption 

Fig. 1 The order parameter profile from the effective 517(3) Polyakov-loop 
theory at the transition point. The correlation length £ can be read off directly. 
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